Self-passivation of vacancies in a-PbO. 
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We introduce a self-passivation of single lead (Pb) and oxygen (0) vacancies in the a-PbO com- 
pound through formation of a Pb-0 vacancy pair. The preferential mechanism for pair formation 
involves initial development of the single Pb vacancy which, by weakening the covalent bonding, 
sets up the crystal lattice for an appearance of the O vacancy. Binding of the Pb and O vacancies 
occurs through the ionization interactions. Since no dangling bonds appear at the Pb-0 pair site, 
this defect has a minor effect on the electronic properties. In such, vacancy self-passivation offers a 
practical way to improve the transport properties in thermally grown PbO layers. 



Polycrystalline Lead Oxide (PbO) is one of the few 
photoconductive materials with a long- more than 60 
years - history of employment in imaging devices. Al- 
though PbO is extensively used as a photoconductor, 
very little is known about the electronic properties and 
charge transport in this material. It is generally believed 
that transport in PbO is controlled by trapping in de- 
fects and that trapping is a cause of low mobility-time 
product [l|. However, the nature of defects is not fully 
understood. Emerging applications of PbO in the direct 
conversion flat panel radiation medical imaging detectors 
revived the interest in studying defects in this material as 
defects can significantly affect imaging performance Q. 

Our recent modeling of the native point defects in a- 
PbO [3] has shown that thermally deposited PbO layers 
should contain a significant concentration of single va- 
cancies due to their moderate formation energies. Single 
vacancies are amphoteric defects appearing in the dif- 
ferent charge states. In the neutral charge state, the 
O vacancy (V°(°))) holds two electrons and forms the 
deep donor level. The neutral Pb vacancy (V Pb (°))) is 
filled with holes, and is a shallow acceptor. It was estab- 
lished that these vacancies prefer to appear doubly ion- 
ized, V Pb ( 2- ) and V°( 2 +), acting as compensating cen- 
tres to each other. Indeed, two compensating vacancies 
have a lower formation energy than neutral ones, such as 
AEf(V Ph ^ AEf(V Ph ^ +y°( 2 +))=0.78 eV 

a. 

The fact that vacancies prefer to appear in charge 
states suggests that we have to consider the ionization 
interactions between them and the formation of a neutral 
vacancy pair V Ph '° instead of two separate compensat- 
ing vacancies. This should further decrease the formation 
energy, approximately speaking, the free energy required 
to insert a defect in a lattice is reduced by the energy 
liberated due to ionization of the donor and acceptor. 
Indeed, the formation of the defect complexes in favor 
over the single defects is often observed during material 
deposition [4]. In this work, we present our study of the 
formation mechanism of V Ph ~° pair in PbO layers and 
its effect on the electronic properties. 



Analysis of the formation of V Ph ~° vacancy pair was 
performed using the density functional theory (DFT) 
available in the Wien2k package [H] which utilizes the full- 
potential augmented plane-wave method. The Perdew- 
Burke-Ernzerhof parameterization [6[ of the generalized 
gradient approximation (GGA) to DFT has been imple- 
mented. We assigned only 5p, 5d, 6s and 6p electrons of 
the Pb atom and 2s and 2p electrons of the O atom to the 
valence states, while the lower energy electrons were in- 
cluded into the core shells. The supercell approach was 
applied in which the supercell was taken of 120 atom 
size (5x3x2 array of the primitive cells shown in Fig. [1] 
(a)). For some calculations, the supercell of 108 atom 
size has been used (3x3x3). The supercell has been opti- 
mized with respect to the internal degrees of freedom, for 
which minimization based on the forces has been applied 
[7j. The following parameters have been used: RK max =7 
(product of the atomic sphere radius and the plane- wave 
cut-off in k-space), and standard Monkhorst-Pack mesh 
of size 4x4x2 in k-space for Brillouin zone. Assuming 
that PbO layers are composed primarily of a-PbO single 
crystals [8j, we have considered it as a model compound. 
The a-PbO single crystal is of the tetragonal symmetry 
(129P4/nmm), its crystal structure is layered and the 
layers are held together due to interlayer orbital overlap 
of the Pb:6s 2 lone pairs [9[. The detailed information on 
structure and parameters used in our investigation of the 
vacancy defects can be found elsewhere [3[. 

Upon formation of the vacancy pair, deficiency of two 
electrons at the Pb vacancy is compensated by two elec- 
trons occupying the O vacancy and dangling bonds are 
closed on each other. Since Pb vacancy state is delocal- 
ized beyond the vacancy site (its wavefunction is spread 
over at least five nearest-neighbours), it is not required 
for O and Pb vacancies to be located at the nearest sites 
to form a pair: schematic presentation of a position of 
the interacting vacancies relative to each other is shown 
in Fig. [1] (b). The electronic interactions between va- 
cancies are calculated as function of separation distance. 
Our analysis has revealed that an extra electron or hole 
added to the pair is delocalized implying that V Ph ~° has 
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FIG. 1: Colour on-line, (a) The crystal structure of a-PbO is 
shown for supercell of size 5x3x2. (b) Scheme demonstrates 
a position of the O vacancy (numbers 1,2,3 and 4) relative to 

vacancy pair. 



the Pb vacancy to generate the V Ph '° 



only zero charge state for which the formation energy is 



only 

AE f(v Ph -° 



E tot {D) - E t ot(X) + /i (0) + M(Pb) (1) 



where E to t(D) and E to t(X) are the total energy of the 
system containing the defect and defect-free, respec- 
tively; /i (0 ) = E tot {0) + n* {Q) and /i (Pb) = E tot (Pb) + 
^(Pb)> are ^ ne chemical potentials of removed atoms 
(O and Pb atoms) where \i\ is the energy required to 
move atom from vacuum to the conditions of growth 
0(Pb) +^( ) = A/ff°(PbO) =-2.92 eV [3] is defined 
for the stable compound). 

To account for ionization interactions between vacan- 
cies we have analyzed the non-interacting and interact- 
ing vacancies. For non- interacting case, V° and V Ph 
have been placed into the same supercell, but in differ- 
ent layers. The interactions between such vacancies are 
absent and each vacancy creates its own localized state, 
and y Pb (°). The formation energy under typical 
growth conditions (/i* Pb ) + /i* Q ) =-2.92 eV) is found to 
be AEf(V°^ + yPb(o)) =5 41 eV Tq initiate the inter _ 
action between vacancies, they have been placed in the 
same supercell and the same layer. The formation ener- 
gies of the vacancy pair AE? (V Ph ~°) as function of the 
vacancy separation are presented in Fig. [2] (a) (in the 
following we assume that A£^(V Pb "°) is saturated at 
position 4). 

The total lowering in the formation energy upon va- 
cancy binding is found to be E^ n d=-1.47 eV (vacancies 
are located next to each other) which suggests that V Ph ~° 
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FIG. 2: The formation energy of the vacancy pair as a func- 
tion of the vacancy separation (V° position is referred to 
positions 1,2,3 and 4 in Fig. [1] (b)): (a) formation of the 
yPb-o vacanC y p a i r> (b) binding of V° to previously formed 
yPb (yO^Pb^ binding of yPb tQ already existing V° 
(V Ph ^°). The formation energy of the single charged vacan- 
cies are calculated elsewhere [J. 



is extremely energetically favored. There are two primary 
parts contributing to Eund'- interaction energy Ei nt and 
lattice relaxation E\ at . The E int is the energy change 
induced by the electron exchange between y°(°) and 
y-Pb(o) caus i n g their ionization to y Pb ( 2 ~) and l/ ( 2 +) 
states which is followed by lattice relaxation E\ at . The 
value of Ei n t can be determined through the difference 
between the considered position and the point of satura- 
tion (position 4 is considered here to be the saturation 
point), is equal to 0.55 eV for V° in position 1. E int 
gradually decreases as oxygen vacancy is moved away 
(see Fig. [2] (a)) thus increasing the formation energy of 
a pair. The magnitude of Ei nt =0.55 eV is foreseen to be 
underestimated because the state formed by the V Pb va- 
cancy is largely delocalized and when V° is placed at po- 
sition 4, it already may interact with the periodic copies 
of V Ph . To eliminate this effect, the supercell of larger 
size has to be implemented which is not computationally 
affordable for us. For all positions depicted in Fig.[T](b), 
Ei a t is constant and equals to 0.92 eV as the charge state 
within V Ph is not changed. 

Although vacancy binding lowers the formation en- 
ergy, its value is still high (> 4 eV) because the Pb 
and O vacancies are considered to appear simultaneously. 
The more realistic case appears to be settling as the 
compensating vacancy originates near already existing 
vacancy. We introduce notation V Ph ^° when V Pb is 
formed nearby existing V°, and y°^ Pb for the opposite 
case. The formation energies are: 

AEf(V°^ Ph ) = E tot (D) - E tot (Pb) + M ( o) (2) 
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E tot (D) -£ to t(0) + /i(Pb) (3) 
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where E tot (Pb) and E tot (0) are the total energy of the 
system containing the single vacancy V Pb or V° , respec- 
tively. 

We have plotted the results on the formation energy of 
the vacancy pairs, y°^ Pb and V Ph ^° , in Fig.[2](b) and 
(c), respectively. The different growth conditions have 
been considered: M( Pb )=0.0 eV and /ij^=-2.92 eV are 
defined for Pb-rich/O-poor limit, while fi^ ph) =-2 .92 eV 
and /i* Q) =0.0 eV for Pb-poor/O-rich limit [3]. There is 
a very important trend to be observed in Fig. [2] (b) and 
(c): the values of AE? (V Ph ^°) and A^(V°^ Pb ) are 
significantly lower than the formation energy of the cor- 
responding single vacancies in their doubly ionized state 
[3[ . Taking into account that binding energy is only E^nd 
=-1.47 eV, this trend indicates that the lattice distor- 
tion induced by the single vacancy reduces the energy 
required to generate another vacancy nearby. The single 
vacancy disturbs the lattice periodicity that alters the 
atomic forces between atoms and weakens the covalent 
bonds around the vacancy site. 

Therefore, the Pb vacancy sets up a lattice to accom- 
modate the O vacancy, thus causing a reduction in the 
formation energy AE^(V°^)- A^(V°^ Pb )=4.52 eV 
(V° is in position 1). Negative AE* (V°^ Pb ) means that 
the lattice containing V Pb tends to "suck in" O vacancy 
and formation of y°^ Pb pair becomes spontaneous. We 
expect the Pb vacancies to be completely passivated by O 
vacancies, as even for O-rich limit the formation energy 
Ai^(V°^ Pb ) is close to zero. Similarly, the formation of 
Pb vacancy in close proximity to pre-existing O vacancy 
is energetically more favorable than single double ionized 
yPb(2-) However, the gain in energy is smaller than in 
the previous case: A£^(V Pb ( 2 -)) - A£^(V Pb ^°)=1.96 
eV and because AE? (V Ph ^°) > 0, the lattice distortion 
induced by O vacancy is not sufficient to make a process 
of origin of the Pb vacancy spontaneous. 

To understand an effect of passivation on the transport 
properties, the next step is to examine alteration in the 
electronic properties induced by the neutral vacancy pair. 
We found, since the uncompensated dangling bonds are 
absent at the V Ph ~° defect site, the vacancy pair has a lit- 
tle effect on the band formation as shown in Fig. [3] The 
vacancy pair does not seems to generate any localized 
state within the band gap, but rather induces a slight 
modification into the band behavior. The top valence 
band and its dispersion are almost unaffected by the ap- 
pearance of the vacancy pair and, therefore the hole mi- 
croscopic mobility is predicted to remain the same (it is 
already low as a result of heavy holes ra^=2.44 mo jllj). 
In the conduction band, the vacancy pair induces a shift 
of the lowest band down on the energy axis. This band 
is a product of closing of the dangling bonds of vacancies 
on each other, i.e. its antibonding orbital. Moreover, it 
shows a slightly weaker dispersion at the M* point that is 
anticipated to cause a reduction in the electron mobility. 

To gain insight the band behaviour we have plotted 
the electron density map for bottom of the conduction 
band (Ec) and top of the valence band (Ey) in Fig. 2] 
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FIG. 3: Band diagram for ideal supercell (black solid line) 
and for the same supercell containing the V ph ~° defect (red 
dashed line). 
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FIG. 4: Colour on-line. The electron density map (e/A 3 ) 
plotted with XcrySDen 12] for 3x3x3 supercell: the cutting 
plane is perpendicular to layers and passes through the va- 
cancy site, (a) the conduction band (Ec+0.076 eV), (b) the 
valence band (Ey-0.076 eV). 



(a) and (b), respectively. The O vacancy was created 
next to the Pb vacancy. The electron density map is 
plotted for the narrow energy range (Ec +0.076 eV) and 
(EV-0.076 eV) to catch E c and E v in vicinity of the M* 
and r points, respectively. For the conduction band, the 
electron density shows the pronounced localization in the 
layer containing V Ph ~° and, in particular, on the atoms 
around the defect site: continuity of the electron network 
is disturbed that should affect the electron microscopic 
mobility. For the valence band, an alteration in the elec- 
tron density is confined entirely to the vacancy site where 
density vanishes. Although such hole in the electron den- 
sity can affect the carrier transport locally in the vicinity 
of the defect, beyond that density remains undisturbed, 
suggesting that the hole mobility is not altered. 
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From our study, one can conclude that Pb vacancies 
in PbO layers are passivated by O vacancies to form the 
neutral vacancy pair. It seems reasonable to assume that 
negative effect of the vacancy pairs on the transport prop- 
erties of PbO photoconductive layers would be minimized 
as the V Pb -° pair can not act as trapping center (dan- 
gling bonds are absent and extra electron or hole is delo- 
calized). Moreover, as pair formation involves merging of 
single Pb and O vacancies (i^nnd =-1.47 eV) which other- 
wise are active traps, it is foreseen to suppress the carrier 
trapping. As the concentration of O vacancies is higher 
than that of Pb vacancies, the large fraction of y°(°) 
is left uncompensated to be filled with two electrons |3j- 
The experimental evidence of n-type conductivity in PbO 
layers [13j supports this observation. Our finding offer a 
practical way to improve the transport properties of ther- 
mally deposited PbO layers. The post-growth annealing 



(in vacuum or oxygen atmosphere) would initiate migra- 
tion of the O vacancies towards Pb vacancies with driving 
forces defined by the vacancy interaction E int =0.55 eV 
and facilitates their merging. This is expected to reduce 
an amount of ionized centres in PbO thus improving drift 
mobilities of electrons and holes as well suppressing re- 
combination. 
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